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This project is under the supervision of Ivan I. Mueller, 
Professor of the Department of Geodetic Science at The Ohio State 
University, and is under the technical direction of James P. Murphy, 
Special Programs, Code ES, NASA Headquarters, Washington, D. C. 
The contract is administered by the Office of University Affairs, 
NASA, Washington, D, C. 20546. 
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1. INTRODUCTION 


The basic purpose of this experiment is to compute reduced 
normal equations from the observational data of the South American 
Satellite Deusification (PC-1000) Network obtained from the Defense 
Mapping Agency Aerospace Center, St. Louis. These reduced normal 
equations are to be combined with reduced normal equations of other 
satellite networks of National Geodetic Satellite Program [Mueller et 
al. , 1973] so as to provide station coordinates from a single least 
square adjustment. 

Details of this network, including instrumentation, are given in 
Huber [1971]. 


2. DATA 

2. 1 Terrestrial Data 

Terrestrial data, which include base-lines, heights and sur- 
vey coordinates of stations, provide the necessary relative position 
constraints between "collocated" stations of BC-4 World-net and the 
South-American Densification Net (Figure 1). Surv^ information 
regarding the observation stations is summarized in Table 2.1-1. 
Constraints used in this solution are given in Tables 2.1-2, 2.1-3, 
2.1-4 and 2.1-5 [Mueller, et al. 1973]. Geoidal undulations 
(Table 2.1-5) are computed by using formula and constants as given 
in [Rapp, 1973]. 
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Fig. 1 South American densification net. 
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Table 2.1-1 



S U R 

V £ Y I 

N F 0 

R H 

A T 1 0 N 

0 

F 

0 R S E 

1 V A T I 0 

N S T A 1 

r I 0 N S 




STATION 

j DATUHI 

S U 

R V E Y 

C 

0 0 

R D I N A 

T E S" 

I MSL® 

INSTR, ^ 

INSTR. 

SOURCE 



1 

1 








HEIGHT"* 



NO 

i NAME 

1 CODE 

1 

LATITUDE 

LONGITUDE 

ELL. H(H) 

t (M) 

(M) 

TYPE 

CODE® 



■ — ^ — 

mm mi 

, m 










3406 

1 

i CURACAO 

1 

1 41 

1 

1 12 

5 

26*843 

291 

9 

45.803 

-4.0 

1 6 .83 

1.25 

PC-tOOO 


3407 

.1 TRINIDAD 

1 41 

1 10 

44 

35.844 

298 

23 

25.652 

237.0 

] 254.80 

1.25 

PC-1000 

1 

3413 

1 NATAL 

1 41 

1 - 5 

54 

56.253 

324 

A9 

57.605 

63.0 

1 36.90 

» 

PC-1000 


3414 

1 BRASILIA 

1 41 

1 -15 

51 

35.540 

312 

6 

2.679 

1059.0 

t 1058.25 

1.14 

PC-IOOO 


3431 

1 ASUNCION 

1 41 

1 -25 

la 

56,192 

302 

25 

15.376 

162.0 

1 149.74 

1.65 

PC-IOOO 


3476 

1 PARAMARIBO 

1 41 

1 5 

26 

54.645 

304 

47 

44,226 

8.6 

1 18.27 

1.25 

PC-1000 


3477 

1 BOGOTA 
1 

1 41 
1 

1 4 

1 

49 

2.379 

285 

55 

35.482 

2586.0 

( 2557.90 

1.25 

PC-IOOO 

2 

3478 

1 MANAUS 

1 « 

I - 3 

8 

44.820 

300 

0 

59,620 

« 

1 83.60 

« 

PC-IOOO 


3499 

1 QUITO 

[ 41 

1 - 0 

5 

50.463 

281 

34 

49.212 

2706,4 

1 2661.80 

« 

PC-1000 


6002 

1 BELTSVILLE 

1 29 

( 39 

1 

39,003 

283 

10 

26.942 

45.0 

r 44.30 

1.50 

8C-4 

1 

6008 

! PARAMARIBO 

1 41 

1 5 

26 

55.325 

304 

47 

42.832 

fl.7 

f 1B.3B 

1.49 

BC-4 

1 

6009 

1 QUITO 

1 41 

1 - 0 

5 

50.468 

281 

34 

49.212 

2706,7 

{ 2682.10 

1,50 

BC-4 


6019 

1 VILLA DOLORES 

1 41 

1 -31 

56 

33.954 

294 

53 

41 .342 

621,0 

t 608.18 

1.50 

BC-4 

2 

6067 

1 NATAL 
1 

1 41 

1 - 5 

55 

37.414 

324 

50 

6.200 

66*7 

{ 40.63 

* 

6C-4A 

1 


♦ INSUPFICIENT DATA 

1 DATUM CODE r 

29 — NAO 1927 
41 — SAD 1969 

2 GEODETIC COORDINATES OF THE INSTRUMENTAL REFERENCE POINT (OPTICAL/ELECTRONIC CENTER.ETC.) 
ON THE LOCAL GEODETIC DATUM 

3 MEAN SEA LEVEL HEIGHT OF THE INSTRUMENTAL REFERENCE POINT 

4 HEIGHT OF INSTRUMENTAL REFERENCE POINT A80VE SURVEY MONUMENT 

5 SOURCE CODE : 

1 — (CSC. 1971) 

2 — (CSC. 1972/73) 

3 — (HUBER, 1971) 

NOTE t ZERO IN THE LAST DIGIT MAY INDICATE THAT THE DIGIT IS UNKNOWN. 



Table 2.1-2 


Relative Position Constraints 


Stations 

Relative 

Coordinates (Meters) 

Weights 
(Vo® ) 

Au 

Av 

A w 

3413-6067 

-48,64 

-289.13 

1258.05 

3.00 

3476-6008 

36.31 

22.94 

-20.80 

3.00 

3499-6009 

0.00 

0.00 

1 

0.00 



100.00 


Table 2.1-3 

Station Position Constraints 



Station 

Coordinates (Meters) 


Weights 


Stations 

u 

V 

w 

Pu 

Pv 

Pw 

6002 

1130 764.85 

-4830 831.87 

3994 704.05 

0.2415 

0.3425 

0.2725 

6008 

3623 241.00 

-5214 233.74 

601 536.05 

0.2212 

0.2591 

0.1163 

6009 

1280 834.24 

-6250 955.94 

-10 800.59 

0.0776 

0.0852 

0.0593 

6019 

2280 627.09 

-4914 543.17 

3355 402.77 

0.1779 

0. 1366 

0.0741 

6067 

5186 397.12 

-3653 933.25 

-654 276.92 

0.2315 

0,2160 

0.1464 


Weights (Va®) based as per variances obtained in OSU WN14 Solution [Mueller 
et al. , 1973] 
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Table 2.1-4 
Chord Constraint 


Station-Station 

Chord Distance 
(Meters) - 

a X 10 

6009-6067 

4734 132.44 

1.00 


The chord length and the sigma were computed from the adjusted co- 
ordinates of the stations from OSU WN14 Solution [Mueller et al. , 1973] 

Table 2.1-5 


Geoidal Undulations and Heights Used in the Constraints 


Station 

NREF^ 

HCONSTR^ 

(m) 

^CON- 

STR^ 

(m) 

Na 

NAME 

(m) 

3406 

Curacao 

-29.19 

-41.02 

' 4.0 

3407 

Trinidad 

-38.57 

194.88 

4.0 

3413 

Natal 

-12.03 

-5.87 

6.0 

3414 

Brasilia 

-9.88 

1021.23 

6.0 

3431 

Asuncion 

11.98 

137.72 

6,0 

3476 

Paramaribo 

-28.31 

-34.02 

6.0 

3477 

Bogota 

10,71 

2551.44 

6.0 

3478 

Manaus 

-7.17 

53.63 

6.0 

3499 

Quito 

16.73 

2682.74 

6.0 

6002 

Beltsville 

-36.90 

-6.73 

2.5 

6008 

Paramaribo 

“28.31 

-33.91 

4.0 

6009 

Quito 

16.73 

2683.04 

6,0 

6019 

Villa Dolores 

22.80 

609.43 

6.0 

6067 

Natal 

-12.03 

-2.14 

6.0 


1. From [Rapp, 1973] 

2. HCONSTR = MSL+NREF+ii.N, where AN is a correction term for the dif- 
ferences of position and size of the ellipsoids used [Mueller et al. , 1973] 

3. Used in Computing the Weights of the Height Constraints 
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2.2 


stellite Observational Data 


Data of South America Densification Net was obtained from the 
Defense Mapping Agency/Aerospace Center, St, Louis, Missouri. This 
data, which is a punched card-deck, is used without any modification. 
No major blunders are detected. (See Table 2.2-1) 

3. THEORETICAL BACKGROUND 
3. 1 Normal Equations for Optical Observations 

A set of reduced normal equations of the form 
NX + U == 0 


is obtained from the optical observations. The symmetric coefficient matrix 
is composed of 3 x 3 blocks of the form [Mueller, 1968] 

- SMkS’- (SM,;')-' MkJ + 

3X3 3 J 1 


if,. = -lm,? m.; 

3x3 i 1 

where 


= B., P,f B<; 

Bij = SR3 i-a.) Rg (-90® +6) 


10 0 
0 -cos6 0 
0 0-1 


S = Polar motion matrix 
and 


k,L Denotes particular ground stations 

i Denotes particular simultaneous event 

I Denotes any ground station participating in an event 
S Is the summation over all ground stations involved in event j 

S fe the summation over all events observed by ground station 

^ k and/or 1 

Weight matrix associated with any particular ground station 

3 x3 


Pi 3 Weight matrix of any observed direction 

3x3 


6 



Table 2.2-1 


Summary of Simultaneous Observations by Line (SA Network) 


Line 

Station-Station 

No. of Pairs 


Line 

Station-Station 

No. of Pairs 

6002-6008 

23 


3406-3478 

14 

6002-3406 

14 


3406-3499 

4 

6002-3407 

n 


3407-3431 

16 

6002-3476 

7 


3407-3476 

19 

6002-3477 

7 


3407-3477 

23 

6008-6009 

10 


3407-3478 

9 

6008-6019 

36 


3413-3414 

29 

6008-6067 

14 


3413-3431 

2 

6008-3406 

26 


3414-3431 

22 

6008-3477 

3 


3476-3477 

15 

6008-3478 

6 


3477-3478 

2 

6009-6019 

7 


3477-3499 

5 

6009-3406 

14 




6009-3407 

6 




6009-3476 

6 




6009-3477 

5 




6009-3499 

9 




6019-6067 

35 




6019-3406 

19 




6019-3407 

38 




6019-3431 

4 




601 9-3476 

19 




6019-3477 

6 




6067-3407 

3 




3406-3407 

9 




3406-3413 

25 




3406-3414 

41 




3406-3431 

53 




3406-3476 

20 




3406-3477 

13 














Finally, the vector of constant terms is expressed by 

= SMki" pj - (SM,;V lmJ xf] 

i i 5 

where as usual, the superscript ('*) denotes initial approximate values. 

3.2 ^^ Constraints*' Contributions to the Normal Equations 
3.21 General 

Since the coefficient matrix of normal equations is singular, a unique 
least squares solution is not possible. A minimal set of constraints to the 
normal equations provides a unique solution [Blaha, 1971]. 

Two alternative definitions exist for the term '’constraints". The 
absolute constraints represent certain conditions which have to be fulfilled 
exactly and with no uncertainties while the relative constraints (or weighted 
constraints) have the same characteristics as the observations. 

In general the contribution of the functional constraint equations 
G(X, U) = 0 

to the normal equations can be found by bordering the normal equations matrix 


N 

C' 


X 

+ 

U 

_ ...u ^ SG 

C 

-P? 




W" 

— u, wiiere k , — 

dX ' 


After elimination bf Kc , it is easy to find 

^ + C'Pc C]X U - C*Pc W" - 0 

or 

[N + N^]X + U + = 0 (1) 

where and are the contributions to the coefficient matrix and constant 
vector of the normal equations due to the application of constraints. The 
quantities N and U correspond to the original normal equations without con- 
straints (Section 3.1). 

After the constraints are added the normal equations will take the 
usual form 

NX + TJ - 0 

and we are in the position to obtain the contribution from a new set of 
constraints. 
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Constraints can be applied between two station k and 1 or to a single 
station. The contribution of these constraints to the matrix N (3 x 3) blocks) 
and U (3 X 1 blocks) can be schematically expressed in two different ways. 

Assuming that the matrix Pc is always diagonal it is possible to 
express: 

(a) Contribution to the normals due to the constraint applied to station k 



(b) Contribution to the normals due to the constraint between stations k and 1 
Pc Ck Cl W** 



These blocks obtained as indicated above for the corresponding case 
will be the only ones computed and added to the original normal equations as 
expressed by formula (1). 
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3.22 Relative Position Constraints 

Relative position constraints are used in order to combine the normal 
equations obtained from various satellite nets and to constrain ’’double" sta- 
tions or closely situated stations of the same net. The expression for the 
combination of normals can be written as follows: 

[N + r^] X + u + u'^ =0 

R R j 

where N and U computed from (2a), (2b), are the contribution to the ori- 
ginal normal equations (NX + U - 0). 

If the relative position (^u° , Av° , Aw° ) of two stations is known, 
along with the standard deviation of these relative positions, the constraints 
can be formed. In this case the functional constraint equations are 

uj - u? = Au*^ 

v° - V? = Av° 

w° - wf = Aw° 

Therefore 

Cj = I ; 

3X3 3X3 3X3 

and 

US = 0 ; u; 

3X1 

where 

1 

0 ^ 0 

0 0 ^ 

^Aw° 



= -I 

3X3 

= 0 because W” = (X^, L°) = 0 
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and 


nJk = I I = p„ 

3x3 3x3 

n 1 1 = I p, I = B, 

3x3 3x3 

1^1 = n"i. = I P (-1) = -p 

3x3 3x3 3x3 

Thus, the diagonal elements of are added to each element of the diagonal 
of the blocks kk and ii of the matrix of the combined normals N, and sub- 
tracted from the diagonal elements of the blocks k i and ik of N. There is 
no contribution to the vector U. 


3»23 Length (Chord) Constraints 


Chord constraints are introduced when scalar information is avail- 
able between ground stations (e,g. , distances determined through high pre- 
cision geodimeter traversing). The functional constraint equation in this case 
is 

(X, Lc ) = 0 


or 


and 


[(U( - Ui)*" + (Vk - Vi)^ + (wk - Wi) ] = Lk 


(3) 


d 

1x3 


cfi 

1 x3 


tk - UT , Vk - VT , wg - Wi 

I lii 


Lfi lii 


u< - u? , vk - V? , w? - w? 


and 


W = 


<yo 


= variance of the chord 


a priori variance of unit weight 


Then the contribution to the normals are obtained by applyii^ (2a) and (2b) 
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as under: 


(Ck)' Pc Ck 


NJk = 

3x3 

1^11 = (cfi)' ft di 

3x3 

= (Ck)' Pc cS 

3x3 

= (Ck )' Pc 

3x3 

= (C°)' ft vf 

3x3 

The first three expressions in the above are added respectively to 
the blocks Nkk: , Ni i and N<i of N; the last two expressions are added 
respectively to the constant vectors Uk and Ui of U. 

3.24 Station Position Constraints 

Station position constraint is used for the purpose of defining the 

origin of the coordinate system. If the station coordinates (u«®, Vk , wk) of 

station k are to be constrained and if the computed (known) variances of its 

known coordinates are then the equations given in section 

Uk Vk Uk 

3,22 are valid by merely deleting the terms with index i, i.e., Au® = ni. 
Av® = Vie ® , AwJ = Wfc . Hence 

N^k = r Ps I = Ps 

3X3 3X3 

where 
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3.25 Height Constraints 

If the geodetic (ellipsoidal) heights of the station k is to be con- 
strained, then 

nJk = (C?)' Ph C? 

3x3 

where 

C? = [coscOk cosXk , cos<Pk sinXK , sin<pK] 

1x3 

and 


1 



Here and X® are the approximate geodetic coordinates and is the 
variance of the height for station ic . 

The constant vector U? can be computed from 

u? = (cHy Ph w“ 

where 

w” = Hk - Hk®, hJ being the approximate height. 

3.26 Inner Constraints (Free Adjustment) 

Even though the selection of a coordinate system is arbitrary in the 
case of a minimum constraint adjustment, e.g. , in the case of ranging, the 
selection of the six coordinates (at more than two stations) to be constrained 
is very critical, since one set of constraints would give a different solution 
than another set. The ’Tsest" solution is arrived at in a coordinate system 
defined through the use of a set of constraint equations called "inner" con- 
straints [Rinner et al. , 1967]. In this sense, ’Tt)est” means resulting in 
the smallest covariance matrix for the unknowns. Covariance matrices may 
be compared by means of their traces, and the inner constraint equations 
are characterized by the property that the trace of the covariance matrix 
obtained with their use is a minimum among those obtained by adjusting a 
given set of observations augmented by a minimal set of constraint equa- 
tions. This property also implies that the mean square uncertainty of the 
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unknowns is smaller when the inner adjustment equations are used. The re- 
sultii^ adjustment is called a "free” one. The functional inner constraints 
equations can be written as 

C^X = 0 

where X is the set of corrections of the approximate coordinates of the un- 
known points and in the most general application when the "best" origin, 
orientation and scale are sought 


1 

n 

HM 


~0 ~ 

I 

3 xa 

3*“ 

WJ. 

-vT- 

0 

r 

3 X3 

D — 

W2 

t ~ 

1 ' * ■ * 

-V2 1 

1 

ci 

- 

0 

-Wj 

0 

U? 

0 

-Wg 

0 

1 

0 ..... 

U. 1 

» 


V? 

-Ui 

0 

0 

vs 

-ul 

1 

0 1 



Vl 

W1 

U2 

V2 

V? 1 



— 





1 • * • • 


The symbols (u? , v? , w? ) denote the approximate coordinates of the 1 th 
unknown point where both the ground points and the satellite positions are 
considered. 

It is also possible to design a set of constraints that will result in 
the "best" solution for only a subset of the points. In the adjustment 
reported here we were only interested in the gound station unknowns Imply- 
that the trace of only that portion of the covariance matrix correspond- 
ing to the ground station unknowns should be minimized, while the variances 
of the satellite position unknowns should not be included in the minimum sum. 
The constraint equations that will produce such a solution have the same 
form as those producii^ the "best" solution for all the points; however, 

3x3 blocks of zeros are inserted into those positions of I which corres- 
pond to unknowns whose variances are not to be included in the minimum sum. 

The inner adjustment constraint equations can be given a geometrical 
interpretation that appeals to intuition. Let X^* denote the set of approxi- 
mate coordinates of the ith unknown point, dXj denote the corrections to 
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these coordinates, and denote the adjusted coordinates, i.e., 

Xi = x; + dXi 

The first set of constraint equations, C^X = 0, is then equivalent to the 
set of conditions 

L dXi = 0 

1 

The geometrical interpretation of these conditions is that the center of 

gravity of all the points mil not change after adjustment, i.e,, 

I; Xj X dXi = 0 
1 

If the center of the system remains fixed, then the cross products Xj x dXi 
reflect rotations of the points around the fixed center. These constraint 
equations insure that the sums of the rotations around all three coordinate 
axes are zero. The corresponding geometrical interpretation is that the 
mean orientation of the system of points will not change after adjustment 
either. 

Thus, the respective equations cjx = 0 and C^X - 0 effectively 
specify the origin and the orientation of the adjustment coordinate system. 

A third "inner adjustment" equation C 3 X = 0 specifies the scale of the 
system. However, this scale equation is only used when the observations 
themselves do not determine the scale. 

A more complete description of the inner adjustment is described 
in [Blaha, 1971]. 

In summary, if the normal eqxiations with the contribution of all the 
constraints (except inner constraints) are represented by 

[N+]^+l^+N®+N”]X + U + tf+Xf+U®+tf=0 (5) 

or 

NX 4 U = 0 

then the inner adjustment can be obtained by bordering the coefficient matrix 
N of the normal equations as 


N (c^y 


_ 

X 


-U 

I 





C 0 




0 
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Upon the addition of any kind of constraint to the normal equations, 
it becomes necessary to consider also its contribution to L V’PV. The de- 
grees of freedom change as well. In order to compute the proper variance 
of unit weight the latter must be taken into consideration. 

4, THE SOLUTION 

Using the specified constraint of Section 3.2, SA-10 solution is com- 
puted. A general information of this solution is given in Table 4-1. The 
solution is obtained by using the general OSUGOP program (Reilly et al. , 
1972]. 

The station coordinates of SA-10 solution are given in Table 4-2 with 
their standard deviations. 

Table 4-1 


General Information on the SA-10 Geometric Adjustment 


No. of observing stations 

14 

Ca (a priori) 

1.0 

No. of degrees of freedom 

649 

EV'PV 

721.76 

O'a (a posteriori) 

2.50 

Number of constraints used: 


Relative position constraints 

3 

Length (Chord) constraint 

1 

Station position constraints 

5 

Height constraints 

14 

Inner constraint defines the 


origin of the coordinate system. 
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Table 4-2 


Cartesian and Geodetic Coordinates 
(Solution SA 10) 



u, V, w Cartesian coordinates in meters (Orientation: u = the Greenwich 
meridian as defined by the B.I. H. ; v = X - 90° (E); w = Con- 
ventional International Origin). 

<P Geodetic latitude and longitude in ai^ular units (degrees, min- 

utes and seconds of arc) computed from the Cartesian coordin- 
ates and referred to a rotational ellipsoid of a = 6378155.00m 
and b = 6356769,70 m. 

H Geodetic (ellipsoidal) height in meters referred to the same 

ellipsoid. 

ot,Qv,crv Standard deviations of the Cartesian coordinates in meters. 

Standard deviations of the geodetic coordinates in seconds of 
arc. 

(Th Standard deviations of the geodetic height in meters. 

a*,Ae,r* Altitude (elevation ar^le), azimuth and magnitude of the major 

semi axis of the error ellipsoid, respectively. Armies in de- 
grees, magnitude in meters. Altitude is positive above the 
horizon. Azimuth is positive east reckoned from the north. 

at,Ab,rb Same as above for the mean axis of the error ellipsoid. 

ac,Ac,rc Same as above for the minor axis of the error ellipsoid. 
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Table 4.2 (Con’t) 


3406 

2251769 . OB 

5.57 
0.24 , 

- 5816902.96 
291 9 43.15 

3.94 

1327210,38 

1 

7.37 

.... 

1 ? 5 26.57 

0,26 

- 46.00 

3.53 



8,88 

“ 32.18 

i 

8.16 ' 





- 5.04 

57.03 

4.83 





- 79,77 

-62 .24 

3.33 



3407 

2979660.94 

7.18 

“ 5513533,18 

4.49 

1181 J 38.43 

6.20 


10 44 35*20 

0,20 

298 23 23.08 

0.25 

165.59 

3.71 



- 3.40 

118.77 

0.31 





20.28 

30.03 

5.44 





-69 ,41 

19.66 

3.37 



3413 

5186345.51 

2.00 

- 3654218.56 

2.10 

- 653016,68 

2.71 


- 5 54 57.55 

0,09 

324 49 55.45 

0,07 

- 4,76 

1.97 



- 3.79 

3.66 

2.71 





- 8.55 

94,23 

2.14 





80.64 

69.9 5 

1.96 



3414 

4114977.26 

. 7.50 

- 4554124.02 

5.71 

- 1732149,84 

7.17 


-15 51 37.37 0.24 

312 6 0.26 

0.27 

1002,05 

4 .88 



2.86 

50,36 

9,10 





- 11.13 

“ 39.07 

5.84 





- 70.50 

126.06 

4.82 



3431 

3093028.20 7.17 

- 4870063.97 

6,08 

- 2710839.99 

10,51 


-25 18 58.25 0.36 

302 25 12.20 

0.26 

130.54 

4.94 



- 5.20 

11,62 

11,19 





-1 .83 

101.79 

7.07 





- 84.48 

- 148.88 

4,05 



3476 

3623275.84 2.14 

- 5214208 .28 

1 . B 7 

601517*56 

2-84 


5 26 52.81 0.09 

304 47 41.51 

0,07 

- 39.47 

1.81 



3.15 

- 6.32 

2.85 





10.83 

84.28 

2.20 





- 78.72 

67.69 

1,79 
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Table 4,2 ( Con ^ t ) 


3477 

1744632.76 

9.75 

- 6114278.13 

6.57 

532213.19 

9,76 


4 49 0.44 

0,32 

285 55 31.56 

0.34 

2542.42 

5.37 



- 2,76 

47.57 

13.26 





“ 43.15 

140.16 

5.79 





46.72 

134.64 

4.90 



3478 

3185743.08 

11 . 19 

- 5514590.21 

9,39 

- 347640,88 

20.43 


- 3 8 43.73 

0.67 

300 0 53.09 

0.43 

36;97 

5.75 



- 0.19 

- 30,23 

23.24 





- 21.96 ' 

59,85 

7.91 





68.04 

59,31 

5,31 



3499 

1280825.47 

2.99 

- 6250950 . 21 . 

2.60 

- 10793.24 

4,01 


~ 0 5 51.25 

0.13 

281 34 46.84 

0,10 

2676.43 

2.61 



2.51 

14.32 

4.0 7 





10.43 

104.78 

2.90 





- 79.26 

90.96 

2.60 



6002 

1130763,54 

2.10 

- 4830831.02 

1.63 

3994704.49 

1.86 


39 1 39.39 

0.06 

283 10 27.00 

0,09 

- 7.30 

1.54 



■ 1.79 

105.75 

2.11 





11.74 

, 15.36 

1.93 





- 78.12 

24.27 

1.52 



6008 

3623239.59 

2.07 

- 5214231.39 

1.79 

601538.31 

2.79 


5 26 53.49 

0.09 

304 47 40,11 

0.07 

- 39.20 

1.74 



3.25 

- 5.65 

2.80 





11.07 

84,98 

2,12 





- 78.4 6 

68.21 

1.71 



6009 

1280825,47 

2.99 

- 6250950.22 

2.59 

- 10793.24 

4,01 


- 0 5 51.25 

0.13 

281 34 46.64 

0,10 

2676.43 

2.61 



2,51 

14,32 

4,07 





10.48 

104.78 

2.90 





- 79.22 

91.01 

2.59 



6019 

2280626.58 

2.43 

- 4914540.71 

2.64 

- 3355401.90 

3,65 


-31 56 34.95 

0.11 

294 53 38.34 

0.10 

603.72 

2.77 



- 24.10 

- 0,19 

3.67 





- 46.19 

117.61 

2*63 





33.99 

72.26 

2.41 
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Table 4.2 (Con*t) 


6067 5186394 . 1 . 9 ? - 3653929.59 2.03 - 654276.72 2.64 

- 5 55 36.71 0.09 324 50 4*04 0*07 - 0.90 1.90 

- 3.90 3.53 2.65 

- 11.21 94.31 2.06 

78.11 74.64 1.89 

In the tables the rotations 44 i{>, and c are about the w, v, and u 
axes, respectively. The scale factor A is in units of ppm. The unit in the 
variance-covariance matrix, for the elements corresponding to the rotations, 
is radian squared. 

5. COMPARISON WITH OTHER SOLUTIONS 

Table 5.1-1 summarizes the transformation parameters between 
SA 10 and South American 1969 Datum> and Table 5.1-2 between SA 10 
and WN14 Solution [Mueller et al. , 1973]. 
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Table 5.1-1 


SAIO -TO- SAD-69 


SCALE FACTOR AND ROTATION PARAMETERS CONSTRAIVEO 


SOLUTIDM FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PARAMTTfRS 



(USING variances ONLY) 


DU DV DW DELTA OMEGA PSI FP? ILON 

METERS METERS METERS (XI .0+6) SECONDS SECONDS SECONDS 

5^.65 16.5;? 36.33 5.80 -0.79 0. 12 0.21 

VARIANCE - COVARIANCE MATRIX 

<yl = 1.13 

0.857D+0? 0.5840 + 01 0.106D+02 -0.447D-05 0.967D-05 -0.2770-05 -0.1260-0^‘ 

0.5B4D+01 0.749D+02 0. 1860+01 0.738D-05 0.564D-05 -O.358D-06 0.1300-05 

0.106D+02 0.1R6D+01 0. 9230+02 0.137D-05 0. 161D-05 -0.644D-05 -0.979D-05 

-0.447D-05 0.738D-05 0.137D-05 0.195D-11 0.721D-14 0.31*50-14 O.llOn-13 

0.967D-05 0. 5640-05 0.161D-05 0.721D-14 0.1B7D-11 -0.266D-12 -0.152D-12 

-0.277D-05 -0.5580-06 -0. 6440-05 0.319D-14 -0.2680-12 0.1360-11 0.440D-12 

-0.1260-05 0.1300-05 -0. 9790-05 O.llOD-13 -0.1520-12 0.440D-12 0.1680-11 

COEFFICIENTS OF CORRELATION 

O.lOOD+01 0.7300-01 0. 1190+00 -0.4010+00 0.764D+D0 -0.2570+00 -0. 1050+00 

0.7300-01 0.1000+01 0. 2240-01 0.7080+00 0. 4760+00 -0.3550-01 0.116D+00 

0.1190+00 0.2240-01 O.lOOD+01 0.I18D+00 0.1220+00 -0,5760+00 -0. 7880+00 

-0.401D+00 0.708D+00 0.1180+00 O.lOOD+01 0.4380-02 0.2270-02 0.7050-02 

0.7640+00 0.476D+00 0. 1220+00 0. 4380-02 0. 1000+01 -0.1680 + 00 -0.8560-01 

-0.2570+00 -0.3550-01 -0.5760+00 0.2270-02 -0. 1680+00 0,1000+01 0.2920+00 

-0. 1050+00 0.1160+00 -0.788D+00 0.7030-02 -0.8560-01 0. 2920+00 0. 1000+01 
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Table 5.1-1 (Con't) 


3414 

3431 

3477 

6008 

6009 

6019 

6067 


RESIDU4LS V 


VI ( 

SAIO 

) 


V2( SAD 

-69 ) 


VI 

- V2 


0.4 

-2.8 

1.3 

3414 

-0.7 

.8.5 

-2.6 

1,2 

-11.3 

3.9 

3.0 

-0.1 

7.8 

3431 

-5.9 

0.3 

-7.1 

9.0 

-0.4 

,16.9 

17.7 

0.9 

B.3 

3477 

-18.7 

-2.1 

-8.7 

36.4 

3.1 

16,9 

-0.2 

0.3 

1.0 

6008 

3.8 

-8.2 

-12.5 

-4.0 

8.4 

13.5 

-0,9 

-0 .3 

-2.6 

6009 

10,0 

4.9 

16.4 

-10.9 

-5.2 

-19.0 

-0.4 

0.5 

-0,9 

6019 

7.4 

-7,3 

7.1 

-7.8 

7.8 

-6.1 

-0,2 

-0.2 

-0.5 

6067 

4.1 

3,9 

7.3 

-4.3 

-4.0 

-7.8 
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Table 5.1-2 


SAIO -TO- WNK 

^ » 3fi: * ){e 4 : « ^ 4 # * 3 ^ 3^ # 


SCALE FACTOR AND ROTATION PARAKFTEPS CONSTRAINED 


SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION PA’JAME-TERS 


(USING VARIANCES ONLY) 


DU 

DV 

OW , DELTA 

OMEGA 

PSI 

EPSILON 

METERS 

METERS 

METERS (Xl.D+6) 

SECONDS 

SECONDS 

SECONDS 

1.63 

-7.78 

-4.92 -0.44 

-0.13 

0.07 

o.io 




VARIANCE 

- COVARIANCE MATRIX 




= 0.92 






0.365D+01 

0.52ID+00 

0.8800+00 

-0. 1710-06 

0.403D-06 

-0. 3390-07 

-0. 1340-06 

0.5210+00 

0.315D+01 

0,4620+00 

0.2440-06 

0,2860-06 

-0.304D-07 

-0.120D-06 

0.880D+00 

0. 4620+00 

0.363D+01 

-0.444D-07 

0.162D-06 

-0.168D-06 

-n. 2770-06 

-0,1710-06 

0. 2440-06 

-0.4440-07 

0.5050-13 

-0. 1550-15 

0,3900-15 

0.3510-15 

0.403D-06 

0.286D-06 

0.1620-06 

-0. 1550-15 

0.856D-13 

-0,8840-lA 

-0.204D-1? 

-0.339D-07 

-0.304D-07 

-0.1660-06 

0. 3900-15 

-0.884D-14 

0.4310-13 

0.740D-14 

-0.134D-06 

-0. 1200-06 

-0. 2770-06 

0,3510-15 

-0.284D-13 

0.740D-14 

0.5350-13 



COEFFICIENTS OF CORRELATION 



O.lOOD+01 

0.154D+00 

0.2420+00 

-0.3980+00 

0.7200+00 

-0.8550-01 

-0.302D+00 

0.154D+00 

0. 1000+01 

0, 1370+00 

0.612D+00 

0. 5510+00 

-0.6250-01 

-0 .2920+00 

0.242D+00 

0.1370+00 

O.lOOD+01 

-0.104D+00 

0.291D+OD 

-0.424D+00 

-0.629D+00 

-0.398D+00 

0,6120+00 

-0.104D+00 

O.lOOD+01 

-0. 2360-02 

0.P35D-02 

0. 6750-02 

0.720D+00 

0.5510+00 

0.2910+00 

-0.2360-0,2 

O.lOOD+01 

-0.1450+00 

-0 .4200+00 

-0. 8550-01 

-0.6250-01 

-0.4240+00 

0.8350-02 

-0.1450+00 

0. 1000 + 01 

0 .1540+00 

-0.302D+00 

-0.2920+00 

-0.629D+00 

0.6750-02 

-0.4200+00 

0.1540+00 

0. 1000+01 
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Table 5.1-2 (Con»t) 


RESIDUALS V 



VU 

SAIO 

) 


V 2 ( WN 14 } 


VI 

- V 2 


3 A 06 

6.6 

-5 .3 

- 14.5 

3406 

- 1.2 

1.5 

3.0 

7.8 

- 6.7 

- 17.6 

3407 

4.9 

3 .5 

- 4.3 

3407 

- 2.1 

- 1.9 

3.1 

7.0 

5.4 

-“^4 

3413 

0.1 

- 0.2 

1.0 

3413 

- 0.1 

0.3 

- 0.9 

0.2 

- 0.5 

1.9 

3414 

-1 .4 

-6 .6 

- 1.2 

3414 

1.5 

7.6 

1.2 

- 2.9 

- 14.2 

- 2.4 

3431 

6.3 

- 5.8 

9. 1 

3431 

- 7.1 

6.7 

- 9.7 

13.3 

- 12.5 

le.p 

3476 

“ 0.8 

0.3 

- 0.3 

3476 

0.8 

- 0.4 

0.3 

- 1.7 

0.7 

- 0.6 

3477 

6.6 

-2 .3 

- 1.5 

3477 

- 7.2 

2.4 

1 .4 

13.8 

- 4.7 

- 2.9 

3478 

8.0 

2 .4 

- 15.4 

3478 

- 22.4 

- 5.8 

45.4 

30.4 

8.2 

- 60.7 

3499 

2.0 

- 0.5 

- 2.8 

3499 

- 2*9 

0.9 

2.9 

4.9 

- 1.4 

- 5.7 

6002 

- 0.5 

0 .7 

1.1 

6002 

0.5 

- 0.8 

-1 .1 

-1 .0 

1.5 

2.2 

6008 

- 0.8 

0.4 

- 0.3 

6008 

0.9 

- 0*5 

0.3 

- 1.7 

0.8 

- 0.6 

6019 

-1 .8 

1 .6 

0.5 

6019 

1.7 

- 1.7 

- 0.5 

- 3,5 

3.4 

1 . 1 

6067 

0.1 

- 0.2 

0.9 

6067 

- 0.1 

0.2 

- 0.9 

0,1 

- 0.3 

1.9 
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